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Abstract. We use a variational mean-field like approach, the coupled cluster method 
(CCM) and exact diagonalization to investigate the ground-state order-disorder 
transition for the square lattice spin-half XXZ model with two different nearest- 
neighbor couplings J and J'. Increasing J' > J the model shows in the isotropic 
Heisenbcrg limit a second-order transition from semi-classical Neel order to a quantum 
paramagnetic phase with enhanced local dimer correlations on the J' bonds at about 
J' c ~ 2.5 . . . 3 J. This transition is driven by the quantum competition between J' and 
J. Increasing the anisotropy parameter A > 1 we diminish the quantum fluctuations 
and thus the degree of competition. As a result the transition point J' c is shifted to 
larger values. We find indications for a linear increase of J' c with A, i.e. the transition 
disappears in the Ising limit A — * oo. 



1. Introduction 

The study of zero-temperature phase transitions driven by quantum fluctuations has 
been a subject of great interest to physicists during the last decade, see [TJ |2J EH H] and 
references therein. For order- disorder quantum phase transitions we basically need the 
interplay between the interparticle interactions and fluctuations. A canonical model to 
study quantum phase transitions is the spin-half Heisenberg model with competing 
interactions in two dimensions. Competition between bonds appears in frustrated 
systems. Besides frustration, there is another mechanism weakening the ground-state 
Neel order in Heisenberg antiferromagnets, namely the competition of non-equivalent 
nearest-neighbor (NN) bonds leading to the formation of local singlets of two (or even 
four) coupled spins. By contrast to frustration, which yields competition in quantum as 
well as in classical systems, this type of competition is present only in quantum systems. 

Recent experiments on SrCu 2 (B0 3 )2 E] and on CaX^OgO IH| demonstrate the 
existence of gapped quantum paramagnetic ground states in (quasi-)two-dimensional 
Heisenberg systems and have stimulated various theoretical studies of quantum spin 
lattices with competing interactions. A famous example for competition in a frustrated 
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Heisenberg antiferromagnet is the spin-half J\ — J 2 model on the square lattice, where 
the frustrating J 2 bonds plus quantum fluctuations lead to a second-order transition 
from Neel ordering to a disordered quantum spin liquid, see e.g. [EU EH HH E21 UHj • On 
the other hand, it has been predicted that frustration may lead to a first-order transition 
in quantum spin systems in contrast to a second-order transition in the corresponding 
classical model EH E3 El • An example for competition without frustration is the 
'melting' of semi-classical Neel order by local singlet formation in Heisenberg systems 
with two non-equivalent nearest-neighbor bonds like the bilayer antiferromagnet [T£|ll9j. 
the J — J' antiferromagnet on the square lattice [T71 120J EI] and on the depleted square 
(CaVO) lattice [221 ESJ- 

In the above mentioned papers [13 HHJ HHJ EDI 1221 EH! the strength of quantum 
fluctuations is tuned by variation of the exchange bonds. Alternatively, the strength of 
quantum fluctuations can be tuned by the anisotropy A in an XXZ model. 

In this paper we study the influence of the Ising anisotropy on the zero-temperature 
magnetic order-disorder transition for the J — J' spin-half XXZ antiferromagnet on the 
square lattice. We use a variational mean-field like approach (MFA), the coupled cluster 
method (CCM) and exact diagonalization (ED). We mention that the CCM, being one of 
the most powerful methods of quantum many-body theory, has previously been applied 
to various quantum spin systems with much success [121 H3 EH [23 EHI EH 123 EE] • 

2. Model 

We consider an anisotropic spin-| Heisenberg (XXZ) model on the square lattice with 
two kinds of nearest-neighbor bonds J and J 7 , as shown in FigCJ 

H = J J2 + s K + As * s J) + Jl E + W + • (*) 

<ij>i <ij>2 

The sums over < ij >j and < ij >2 run over the two kinds of nearest-neighbor bonds, 
respectively (cf. FigHJ). Each square- lattice plaquette consists of three J bonds and one 
J' bond. We consider antiferromagnetic bonds J' > J > 0, i.e. there is no frustration 
in the model. In what follows we set J = 1 and consider the Ising anisotropy A > 1 
and J' as the parameters of the model. Since there is no frustration the classical ground 
state is the two-sublattice Neel state. 

3. Methods 

3.1. Variational mean- field approach (MFA) 

For the square-lattice antiferromagnet (J = J') the ground state is Neel ordered. 
The corresponding uncorrelated mean-field state is the Neel state \(pMFi) = | T}| J. 
)| T)l I) • • • • I n the limit J' — > 00 and for finite A the ground state approaches a 
rotationally invariant product state of local pair singlets (valence-bond state) \4>mf 2 ) = 
IlieAil ti)l U+i) ~\ U)\ h+x)} /V2 , where % and i + x correspond to those sites 
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Figure 1. Illustration of arrangement of bonds in the J — J' model on the square 
lattice (Eq.JTJ): J - dashed lines; J' - solid lines; A and B characterize the two 
sublattices of the classical Neel ground state. 



which cover the J 7 bonds. In order to describe the transition between both states, we 
consider an uncorrected product state of the form[T7J HH] 

l* mr ) = jj— L=[| Ui+&) -t\ uh+i)}- (2) 

The trial function \^ var ) depends on the variational parameter t and interpolates 
between the valence-bond state \4>mf 2 ) realized for t — 1 and the Neel state |0mfi) 
realized for t = 0. By minimizing E var = (ty var \H\ty var ) with respect to the variational 
parameter t we obtain 

En* = f -2iA (J' 2 + 3 A 2 J' + 9A 2 ) for J' < 3A 
AT 1 -§J'(A + 2) for J' > 3A . U 

The relevant order parameter describing the Neel order is the sublattice magnetization 

*, /t , » it v f iA/l-l^'/SA) 2 forJ'<3A 
M.H*~l^l*~> = |jj V forJ ^ 3 A. (4) 

M s vanishes at the critical value = 3AJ. The corresponding critical index is the 
mean-field index 1/2. Eq. J3J) may be rewritten in terms of M s as 

Kar/N = -~ J'A - ±JVl-4M? - -AM 2 . (5) 

We expand -E„ ar up to the fourth order in M s near the critical point and we find a 
Landau-type expression, given by 

Evar/N = -ij'(A + 2) + i(J' - 3A)M S 2 + i J'M S 4 - (6) 

o z z 
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3.2. Coupled cluster method (CCM) 

The CCM formalism is now briefly considered, for further details the interested reader 
is referred to Refs. fZT\ EH] • The starting point for the calculation of the many-body 
ground state is a normalized reference or model state |$). For our model the appropriate 
reference state is the Neel state. To treat each site equivalently we perform a rotation 
of the local axis of the up spins such that all spins in the reference state align in 
the same direction, namely along the negative z axis. After this transformation we 
have |$) = | |)| J,) | |)| I) . . . . Now we define a set of multi-spin creation operators 
Cf = sf , sfsf , sfsfsf , . . . . The choice of the Cf ensures that {<&\Cf = = Cj\4>) , 
where Cj is the Hermitian adjoint of Cf. The CCM parametrizations of the ket and 
bra ground states are then given by 

|v&)=e 5 |$), S = J2SiCf, (7) 

= {§\Se- s , S=1 + J2^i C i- (8) 

The correlation operators S and S contain the correlation coefficients Sj and Si which 
have to be calculated. Using the Schrodinger equation, H\^f) = E\^f), we can now write 
the ground-state energy as E = (<&\e~ s He s \<&) . The order parameter is calculated by 
M s = -(*|sf|*) . 

To find the correlation coefficients Si and Si we have to require that the expectation 
value H = (^f\H\^f) is a minimum with respect to Si and Si. This formalism is exact if 
we take into account all possible multispin configurations in the correlation operators S 
and S. Of course, this is impossible for our model. Thus we have to use approximation 
schemes to truncate the expansion of S and S in the Eqs. (JZJ) and (jHJ). The most common 
scheme is the LSUBn scheme, where we include only n or fewer correlated spins in all 
configurations (or lattice animals in the language of graph theory) which span a range 
of no more than n adjacent lattice sites. 

To improve the results it is useful to extrapolate the 'raw' CCM-LSUBn results 
to the limit n — > oo. Although no exact scaling theory for results of LSUBn 
approximations is available, there are empirical indications 123 EZl EH] of scaling 
laws for the order parameter for antiferromagnetic spin models. In accordance with 
those findings we use M s (n) = M s (oo) + ai(l/n) + a2(l/n) 2 to extrapolate to n — > oo. 
Vanishing M s (oo) determines the critical point J' c . The values for J' c obtained by 
extrapolation of the LSUBn results for M s are, however, found to be slightly too large 
[T7] . We may also consider the inflection points of the M s (J r ) curve for the LSUBn 
approximation, assuming that the true M S (J') curve will have a negative curvature 
up to the critical point. We might expect that (for increasing n) the inflection point 
J' in j approaches the critical point J' c . Thus determining the inflection points for the 
LSUBn approximation again we can extrapolate to the limit n — > oo using a scaling 
law J' in f(n) = J- n j(oo) + 61(1/72) + b 2 (l/n) 2 and interpret J' in f{oo) as the critical value 
J' c . 
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3.3. Exact diagonalization (ED) 

In addition to the variational mean-field approach and the CCM we use ED to calculate 
the order parameter. We consider finite square lattices of N = 8, 10, 16, 18, 20, 26, 32 
sites and employ periodic boundary conditions. The relevant order parameter for 
finite systems is the square of the sublattice magnetization M 2 , here defined as 
= {[jj J2f=i T i s i} 2 ) wrt h the staggered factor r i6 A = +1, r igjB = — 1- For the 
finite-size scaling of M 2 we use the standard three-parameter formula (2H1 C2D1 ED E2] 
M s 2 (iV) = M s 2 (oo) + ciN~ 1 / 2 + c 2 N~ l . The critical value J' c is that point where M 2 (oo) 
vanishes. 

4. Results 

To illustrate the behavior of the oder parameter in dependence on J' we present M S (J') 
calculated by CCM (FigJ2J) and by ED (FigJHJ) and the resulting extrapolated values for 
a particular value of Ising anisotropy A = 2. Notice that corresponding results for the 
isotropic Heisenberg case (A = 1) can be found in Ref . [T7] . 

From the extrapolated order parameters one gets the critical values for A = 2: 
J' c = 6.46 (CCM) and J' c = 4.97 (ED). The mean-field value is J' c = 6. The inflection 
points of the M S (J') curves in Fig|2 are J' inf (n) = 5.57 (LSUB2), 5.42 (LSUB4), 5.26 
(LSUB6), 5.11 (LSUB8) leading to an extrapolated value of J' c = J' in Joo) = 4.66. As 
mentioned above, the extrapolation of the CCM results of the order parameter tends 
to overestimate the critical value and yields the largest J' c . This is connected with the 
change of the curvature in the M s -J' curve in the vicinity of the critical point, cf. Fig. 
I2J Therefore the critical value J' c taken form the inflection points seems to be more 
realistic. Obviously, the difference in J' in ^ between the LSUBn approximations is small 
and the extrapolated value is quite close to the value for LSUB8. This statement holds 
for all values of A. E.g. for A = 4 one finds J' inf = 10.37 (LSUB2), 10.29 (LSUB4), 
10.13 (LSUB6), 9.94 (LSUB8) leading to an extrapolated value of J' c = J' inf {oo) = 9.41. 

Our results for the critical point J' C (A) obtained by MFA, CCM and ED are collected 
in FigEJ We find that the CCM results obtained by the extrapolation of the order 
parameter are in good agreement with the MFA data. On the other hand, there is 
an excellent agreement between the CCM results obtained by the extrapolation of 
the inflection points and the ED results obtained by the extrapolation of the order 
parameter. Clearly we see indications for a linear increase in J' c as predicted by mean- 
field theory. 

We mention that the curves shown in Fig0] cannot be extrapolated to A < 1. 
Similar to the effect of the Ising anisotropy (A > 1) one rather expects an increase of 
J' c due to XY anisotropy, i.e. for < A < 1. Indeed for the pure XY J — J' model 
(A = 0) the critical value was estimated to J' c = 4.56J[33J. 
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Figure 2. Sublattice magnetization M s versus J' for A = 2 using coupled cluster 
method (CCM), sec text. 




Figure 3. Order parameter versus J' for A = 2 using exact diagonalization of finite 
lattices of different sizes N, see text. 
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Figure 4. The critical points J' c versus anisotropy parameter A using mean-field 
approach (MFA), CCM with extrapolation of the order parameter (CCM I), CCM 
with extrapolation of the inflection points (CCM II) and exact diagonalization (ED). 

5. Summary 

We have studied the zero-temperature magnetic ordering in a square-lattice spin-half 
anisotropic Heisenberg (XXZ) model with two kinds of nearest-neighbor exchange bonds 
J and J', see FigJTJ In particular we discuss the influence of the Ising anisotropy A on 
the position of the quantum critical point J' c separating the phase with semi-classical 
Neel order (J 7 < J' c ) and the quantum paramagnetic phase without magnetic long-range 
order (J 7 > J' c ). For this we calculate the order parameter within a variational mean- 
field approach, the coupled cluster method and exact diagonalization of finite lattices up 
to iV = 32 sites. We find in good approximation a linear relation J' C (A) oc aA (A > 1) 
with a ~ 2.3 . . . 3.0. This result can be attributed to the reduction of quantum spin 
fluctuations with increasing Ising anisotropy. In the pure Ising limit (A — > oo) the only 
remaining z — z terms in the Hamiltonian commute with each other, i.e. no quantum 
spin fluctuations are present and, consequently, the critical point disappears in the pure 
Ising limit. 
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